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§1. TS a series of papers”) dating back to 1929 the 
operational calculus, as originally developed by O. 
Heaviside, was extended in different directions. Whereas 
Heaviside applied its methods mainly to problems con- 
nected with electrical circuits with constant parameters 
(independent of the time) which were initially at rest, the 
method could be extended to arbitrary initial conditions. 
But also differential equations with variable coefficients 
could be attacked with this method and the properties of 
the solutions studied. Many new relations between different 
functions such as Bessel functions, spherical harmonics, 
Laguerre’s and Hermite’s polynomials, and others, could 
thus be obtained in an extremely simple way. Also the 
method could be extended to a “simultaneous opera- 
tional calculus ” where two or more variables are treated 
operationally. Again, following Heaviside, the method 
appeared a powerful weapon for evaluating integrals. 
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It is well known that in essence Heaviside’s operational 
calculus and the extensions obtained later on are based 
on the Laplace transform. This point was stressed by 
J. R. Carson) and on several occasions by ourselves ©), 
But it is, perhaps, not so well known that O. Heaviside 
himself plainly and clearly points this out in his book 
‘Electromagnetic Theory,’ iii. p. 236 (London, Benn 
Brothers, 1922). 

The great heuristic value, however, of the symbolic 
calculus does not seem to be generally realized yet, 
although a set of simple rules have been developed which 
make the use of it a relatively simple matter. 

It is the purpose of the present note to apply the sym- 
bolic calculus this time to some problems connected with 
the theory of prime numbers; the ease and simplicity 
with which the results are obtained may, we hope, again 
show the great heuristic power of it. But the use of the 
symbolic calculus in this field shows another great ad- 
vantage over the usual methods in so far that the develop- 
ments can most easily be depicted graphically, leading to 
a set of graphs which, we think, may considerably clarify 
the properties of many discontinuous functions belonging 
to the theory of numbers. This possibility is brought 
about by the fact that the operational “image” of a 
discontinuous “original” is a continuous function. 
As, moreover, the “original” fully determines the 
‘image ” and vice versa), we can deduce the properties 
of the discontinuous “ originals ” from those of their 
continuous “images.” The great gain thus obtained 
will be obvious. (A further general extension of the 
symbolic calculus is in course of publication.) 


§2. If we have given a function h(x) which we shall 
call the “ original ” we obtain its “ image ” f(p) with the 
aid of the integral 


fir)=p\ emhayde,. .  . O 
which relation we shall denote by 
Spek), 2 2 wwe. 2) 


assuming at the same time, as has been usual so far in 
the symbolic calculus, all our originals to be zero for 
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z<0. Thus, quoting some former results), it is easily 
shown with the aid of (1) that 


Pa 
p—a ; er, . . . . . . . $ 3 . . (3 a) 
1 
— =a, o 8 8 ew ee lw ll 3b 
i (3b) 
l a 
p—i =e—], . . . . . . . . . ° (3 c) 
ASA, (A=arbitrary constant) . . (3d) 


1 
be ==logz+C, (C=Euler’s constant) . (3e) 


—log (p—1)==Bi(z) = | as, su ep. a oe Gf 


to which we now add 


log (ap-+-8)==log B+ t ds=log ß—Ei(— Ea), (4) 


which latter relation follows also at once from (1). 

A very important function, of which we shall make 
considerable use further on, is the unit function, which we 
denote by U(x—a) and which is zero for x<a, where a is 
a positive constant, and which equals unity for x>a. 

Its image can, for Re p>0, at once be obtained from (1) 
as follows : 

p| e- rzU(x—a) de=p | e-Pt 1 .dz=e-®, 
0 a 


or, for short: 
e~?=U(a—a). e sw ew ee (5) 


The data from the symbolic calculus here given suffice 
to derive not only the “ Prime number theorem,” but 
also several other relations of the theory of numbers. 


§3. Let us first consider the step function [x], where, 
as usual, the brackets mean the greatest integer contained 
in the real variable x. This function is depicted in fig. 1. 
As indicated by the shading it can be considered as the 
superposition of the unit function U(z—1), jumping at 
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x=1, the unit function U(z—2), jumping at x=2, etc. 
Thus, according to (5), we have 


e-P+e-?4e-%1 |... = [e], 


or 


sH} Re p>0. ... . (6) 


Fig. 1. 


Br ee KW 
Z| 


Fig 
3 


Fererceetey 


It is, further, easy to find the image of the saw-tooth 


function of fig. 2 given by x— [z], for, with the aid of (3 b), 
we have at once 
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or, if we subtract } from both sides (making the average 
zero), 


2 4 1 d 4p 
= al D aa ee, ee 2E 
eof] pi = Eloi) O 
which function of p plays an important rôle in the higher 
theory of the I’ functions. It is here derived, however, 
in an extremely simple and natural way as the image 
of the saw-tooth function of fig. 2. 


$4. Consider next the zeta function of Riemann, C(p), 
as a function of p, and let it be the image of an original 
h,(x) to be obtained. 

Now ¢(p) is defined for Re p>1 by the series 


1 1 1 
Up=itota ts: -3 . . . (9) 
which we can write as 
tlp) =e- rig 14 e—plog?_4 g—P logs + ibre s (9 a) 


But according to (5) e-?!°" is the image of U(x—log n), 
i.e., the function which jumps from zero to the value unity 
at x=logn. Thus the original of ¢(p), see fig. 3, is the 
function which jumps a unit step at x—log 1, x=log 2, 
x=log 3, etc., and therefore we can write down at once 
the very fundamental relation 


ilp) = [e] : Rep>l. . . . (10) 


It is of interest to compare this symbolic relation with 
the one obtained above : 


eee , Rep>0 . . . (6) 


e?— 1 


showing that ¢(p) plays the same rôle with respect to e” 
1 

e— l1 

The fact that (10), which is here derived for the first 

time, could be obtained (like (6)) with the help of the 

symbolic calculus without practically any additional 

calculations again shows the great simplicity of the methods 


plays with respect to x. 


as the function 
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employed. From (10) with the aid of (1) we can now 
write at once the integral 

Up)=p | ele] dx, 
or, replacing e by s, 


ww)=n | EL, eye sD 


Fig. 3. 


1 fi i 2 3 
—— y 


} t f l 
log1 log2 tog? log4 ete. 


defining, like the basic series expansion (9), {(p) for 


Re p>1. 
But again, as in (7), we can from (10) derive a saw- 


tooth function by subtracting from (10) the exponential 
(see 3 a) 
Pp e? 
Per 


1 
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and thus obtain, see fig. 4: 


poi iosi], Rep>l, - . (12) 
which is the analogue of (7). Now it is of interest to note 
that the natural construction of the saw-tooth function 
e’—[e*] in (12), which oscillates between zero and unity, 
just annuls the one and only pole of {(p) at p=1, because 
it is known that the left-hand member of (12) has no 
poles). Again, from (12) with the aid of (1) we can now 
write at once the integral 


oh tol=e | mete} de, Rep>t, 


Fig. 4. 


aa ~ Sp) e% fe] 


Teta 
AV 


0 | 1 | 
log2 log3 log4 tgs etc 


or, replacing e” by s, 


P s-s] 
Zatorej, pa o 2a 
which is already valid for 0< Rep, because (12a) re- 
presents the analytical continuation of the left-hand 
function which now has no poles. 

Thus (12a) extends (11) in a welcome way in so far 
as it converges already for Rep>0. Hence (12 æ) provides 
us with an integral along the real axis which defines {(p) 
also in the critical strip 0<Re p<1, in which all the 
much discussed zeros of the function are situated. 

It is also an easy matter to construct the original of 


1 
(1— g=) t0), 
because it is equal to 


E (IPP Āe o. (18) 


n=1, 2.3... 
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and therefore represents the discontinuous square-topped 
function of fig. 5, which jumps up unity at x=log (2n+1) 
and down unity at x=log (2n), and thus oscillates again 
between 0 and +1. Again the image in (13), representing 
an original which thus remains finite everywhere, remains 
itself finite at p =l. 


§5. In order now to derive in a most simple way the 
“Prime number theory ’’ and some other relations of the 
theory of numbers we note, fig. 3, that 


h(x) = [e*], 


Fig. 5. 


(- 57) 800) 


t 
1 
y2  log3 tog4 1095 ete. 


the original of ¢(p), is everywhere contained between the 
two limits y, and y,, where 


y,(4) =e@—1, 


Y2{t)=e*, 
or operationally 


- . « (14) 


which two functions of p we can consider as a first rough 
approximation to ¢(p), but where full account is taken of 


its one pole at p=1. The first expression, is also 


1 
p—l 
obtained when in (9) we replace the summation by an 
integration as 


A 
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An even better approximation, y(x), lying between y,(z) 
and y(x), is at once evident from a consideration of the 
original h(x) of {(p-+1). For, see fig. 6, t(p+ 1), according 


Fig. 6. 


tg" Eo) 


to (9), is given by 


ESIE 5 ori =le-? og 14 1e—p log? the-P logs cans 
(15) 


the original k(x) of which is represented by the broken 
line of that figure, because, according to the rules of the 
symbolic calculus expounded above, the original of (15) 
is the function which jumps unity at x=log 1, one-half 
at «=log 2, one-third at x=log 3, etc. The ordinate of 
this function at the jump at x=logn is therefore 
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I+ ; n +... 7 It further tends towards a straight 
line. Now the definition of the Euler constant C =0:5772 
is 
; l1 1 1 
Lim Itstet--- Fogn } =C, 

and hence the function h,(x) of fig. 6 tends towards x+C. 
The two corresponding functions y, and y, belonging to 
fig. 6 derived from (14) are operationally 


, : Bee os, 
P 


ee n paar: carl cee 
(p+1)—1 p (p+1)—1 
the originals of which, as shown in fig. 6, are 


y=z=1, 
P 


_ 1 
=a+1=~- +1, 
A P 


The better approximation is therefore y,;=«-+-C, or, in 
the p field, 


1 
=g- = - +C. 
Ys D 
Thus the better approximation in fig. 3 becomes 
a d 
Yz = p—i +C, 
or Yz =e] + C, 


lying, as shown in fig. 3, between y; and ya 


§6. We have thus found two rough approximations 
for (p), viz., 


a : 
Y= —— and are . » . (14) 


p—l 

and a better one, 
a, 

Y3 res hie . . . . . . (16) 


The last expression, moreover, confirms the known 
result (5) ; 


. 1 
re: (t=) =C, 
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so that the approximations used here take full account 
of the only pole of ¢(p). It is of great interest, as we shall 
see further on, that these rough approximations suffice 
already for the derivation of several asymptotic theorems 
of the prime numbers, whereas the irregularities in the 
functions h(x) and h,(x) of fig. 3 and fig. 4 and of the later 


Fig. 7. 


Í Š seplogN y~ 
ho+log (p= eP OIN eT y R 


figures must therefore be produced by the more detailed 
properties, such as the zeros, of (p). 


$7. We are now ready to derive the prime number 
theory. We therefore consider h(x), the original of 
log {(p), for which according to (14) we have the rough 
approximations 


Yı = log (1) and Y= log (574) . (17) 
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and the better one 
1 
= — ow « « = (18 
Yz = log (t0) (18) 


Now the well-known Euler expression relating ¢(p) 
to the ordinary numbers n and the primes N is 


1 1 
= ts. oe, 2 
SC Zw (cI) (19) 


or, taking logarithms, 
log tlp) = > e~Plog NI Je—2p log Nt Le—Sp log Nt AA 


N=primes 
(20) 


The original k(x) of the second member of (20) is again 
easily constructed (see fig. 7); for, taking the first 
prime, viz. 2, first, it means, as depicted at the bottom 
of the figure, the function which jumps unity at z=log 2, 
one-half at x=log 22, one-third at «=log 3%, etc. This 
contributing part due to the first prime 2 is represented 
by the top thin curve at the bottom of fig. 7. This 
accounts for the first prime, viz. 2. We can again con- 
struct a similar function for the second prime, viz. 3, 
for the third, viz. 5, etc., and add all the results. Thus 
the thick curve of fig. 7 is obtained. Now, according to 
(17) and (18) we have for log ¢(p) the approximations 


nog), 


nle) = 1g (2), .. . (21) 
yale) =log( -+ +0), 


which, with the aid of (3), are seen to be equivalent to 
yı(z)=Ei(x), 
Ya(x)=Ei(x)—log x—C, . (22) 
Y3() =log (1—C) — Ei { —«(5 —1)} +-Ei(z). 


These functions are also depicted in fig. 7, which shows 
the very close approximation to the actual discontinuous 
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function, even for small values of x. Moreover, for 
xz l all three functions of (22) tend towards Ei(x). 
But according to the construction of the discontinuous 
curve of fig. 7 the ordinate at z)=log Ny (if we neglect 
for each prime all the jumps following the first one, which 
we shall call the subsidiary jumps) is equal to z(N,), 
the number of primes smaller than No, because at each 
prime the function makes a unit jump. Hence for big 
values of x) we have 


n(No) X Ei(x,) =Ex(log No), 


z -3e s 
and, as Ei(x) 5) T ds =li(e*), 
where * ds 
li(a)= | — 
ile) | o logs’ 
we obtain the asymptotic prime number theorem : 
(NYUN), 2 « = + » = (28) 


or, in words: the number of primes under N is asymptoti- 
cally given by the logarithmic integral of N. 

It may be noted that this asymptotic result is obtained 
from all three approximations (21). Their difference 
becomes only of importance at the lower end of the curve 
of fig. 7, where, as could be expected, y, shows the optimum 
approximation. 


§ 8. In order to derive another theorem of the theory 
of numbers we shall next consider the original h,(x) of 


log {(p+-1). ; 
From (20) we first obtain 


log &(p+1)= Z a 
i 


N= primes S 
Neti 


=E{ eres Xf gh ge a eN agree oF } ; 


N 


and the contribution by each prime to the original can, as 
shown in fig. 8, at once again be read from (24). For, 
taking the prime 2 first, we obtain the upper thin curve 
of the upper drawing of fig. 8. It jumps unity at x=log 2, 
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one-eighth at log 2?, etc., and similarly for the other primes. 
We thus obtain the thick curve of fig. 8, where the top 
part of the curve is reproduced on a bigger scale in the 
inset. 


Fig. 8. 


EPON e POIN $ e3plogn $ 
N7? N? 


hy (xe log Spt) = £ 


Naprimes N 


0 1 2 3 4 5 


— X 


Now from our approximations (17) and (18) to ¢(p) 
we have as approximations for log {(p-++1): 


1 
4,(x) = log -, 
x E 
1 
(£) = 2+1}, 
y (2) = log(5 + ) 


nle) “log (5, +0), 
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the originals of which are easily obtained as : 


y(x) =log +C, 
Y(x)=log x+C—Ei(—z), 


ya(a)=log «-+C—Ei(—&), 


(25) 


which three functions are drawn in fig. 8, which shows them 
to be again very good approximations. 

If we neglect for the moment the subsidiary jumps, 
according to the construction of fig. 8, the value of the 


a i. e., the sum 


of the reciprocals of all primes up to Ny. But it follows 
from (25) that this value must approximately be equal to 
log %+C or log (log No) +C. We have in this way 
neglected the contribution due to the subsidiary jumps. 
This total subsidiary contribution is separately depicted 
in the bottom figure of fig. 8, and it tends towards 0-315... 
as a numerical summation shows. 
Hence we obtain the second prime number theorem : 


N 
function at the abscissa x =log Ny is 3 
2 


No 
X L log (log Ny)+-C—0-315...., 
nN 


or 

: 1 1 1 1 1 l 

Lim S 3+3 +g tg +T y 8 log No) } 
~C—0315..., . . (26) 


where the first member contains the reciprocals of all 
prime numbers up to Ny. Expressed in words : 

The sum of the reciprocals of all primes up to Ny minus 
the logarithm of the logarithm of No asymptotically equals 
C—0-315 . . ., where C is Eulers constant. 

It is of interest to compare (26) with 


Lim d i+s+e+yt es + —beN } xO, 
(27) 
where the reciprocals of all integers occur. Apart from 
the additional constant —0-315... formula (26) is of the 
same form as (27), but log Ny in (27) is replaced by 
log (log Ng) in (26). 
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§9. In order to derive still another theorem we shall 


consider the original of st, As approximations 


tip) 


we obtain at once from (21) 


. 1 
uT 

a 1l 1 
y) = z 2 


« 4d 


the originals of which are 


(28) 


Again, from (19) we obtain 


Ee Ho z” z log N. {e7? le N~} ¢-2plog N 
=primes 
te 8p log Ny ee a}. ee (29) 


The original h,(x) of (29) can again at once be read from 
the form of (29), for the contribution to the total ex- 
pression by each prime is a step function like fig. 1 
the steps of magnitude log N occurring at x=log N i 
æ=2 log N ...,etc. The result is depicted in fig. 9, where 
also the approximations (28) are shown, and again a very 
goodagreementis obtained. All three approximations (28) 
tend, for x>1, towards e?. Neglecting again the sub- 
sidiary steps, which for z>1 is permitted, we thus obtain 
the asymptotic theorem 
x log Nem =ele —N,, 


N=2, 3,5 ..No 
or log (2, 3,5, 7, 11... Nod =No « « (30) 
or in words : 


Asymptotically any prime equals the sum of the logarithms 
of all smaller primes. 
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Fig. 9. 


._ Xip -plogN 4 ~3plogN 
hstx, =- g= logn [er +e POON , 5 3p log +--- 
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The simple construction of (29), which can also be 
written as 


_ tp) 1 
tp) TANEN. para. 


can be shortened with the aid of (6), yielding the opera- 
tional relation 


o'(p) 
Ferme. [ggg]: «0 


§10. Finally we consider the original (x) of — (p+) : 
S(p-+1) 


In the same way as before we obtain as approximations 


J 
neos , 
1 l 
a 
Yx) p pHi 
pazl, 
p pe+l 
the originals of which are 
y(x) =a, 
Yo(2) =2-+e-*—I, 2. a g 7 (32) 


y3(%)=a-+C(e~ © —1). 
Now we easily obtain 


_ E+!) 
(p+1) 


1 1 
S | . > e-plogN ep 2p log N 
= Z logN { xe + xe tenh, 
the original h,(x) of which is shown in fig. 10, which was 
obtained as follows :—The contribution of the first prime, 
viz. 2, consists as depicted again at the bottom of the 
figure in a jump of magnitude 4 log 2 at x=log 2, then 


N = primes 


again a jump of a log 2 at ~=2 log 2, etc., and similarly 


for the other primes. Thus approximately a straight 
line is obtained from primes only, the approximations 
Yı Yz and y; of (32) again being shown. An even better 
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approximation than y, is given by the linear relation 
Ya=x—C, as indicated in the figure. Omitting again the 
subsidiary jumps we thus get the asymptotic relation 


‘slog N 
- N X %=log Ng, 


Fig. 10. 


-p+ _ =plgN 2p 
heb) == pa =2 log N Sta + 


Naprimes 


T 


or 22,38,58,72... NoNo, 


where, again, the product contains all the primes up to 
Nọ. In words: 
Asymptotically any prime equals the productof the Nth roots 
of all smaller primes N. 
3P2 
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In conclusion, the very good fit of all the functions 
considered obtained from rather rough approximations 
of the {(p) function and its derivatives, wherein mainly 
the pole p=1 is considered, and derived from ¢(p)==[e*], 
is noteworthy. Also the fact that the asymptotic 
theorems can thus be derived in an extremely simple way 
with the aid of the symbolic calculus may again be a 
proof of the great heuristic value of the latter. 


Finally, I wish to thank Mr. C. C. J. Addink for con- 
structing and drawing the diagrams. 
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Introduction. 


IDHE behaviour of electrons in chlorine and bromine 
has already been described in various publications t, 
and the present communication gives an account of an 
investigation of the motion of electrons in iodine, thus 
extending our knowledge of the electrical properties 
of the halogens. 
* Communicated by Prof. V. A. Bailey. 
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